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1.  Introduction 

Operational  performance  of  Naval  radars  is  routinely  checked  by 
measurement  of  the  vertical  coverage  diagram  (VCD).  Comparisons  of  returns 
from  a  calibrated  target  with  the  VCD  facilitates  the  detection  of  any 
degradation.  This  may  be  in  the  form  of  a  lower  average  detection  range 
or  'holes'  in  the  vertical  coverage.  The  former  may  result  from 
electronic  degradation  or  transmission  line  losses,  while  the  latter  may 
result  from  antenna  damage  or  multipath  effects  -  these  being  determined 
by  sea  state  and  choice  of  antenna  height  or  operating  frequency. 

In  reference  1,  computer  programs  were  described  which  calculated 
(i)  the  radar  return  from  a  target  flying  a  specified  height/range 
profile  and  (ii)  the  probability  of  paint  for  fluctuating  and  non¬ 
fluctuating  targets.  Refinements  to  the  model,  and  the  theoretical  basis 
of  the  algorithms  were  outlined  in  reference  2. 

Comparisons  have  been  made  between  the  output  of  these  programs 
and  the  measured  returns  in  the  RAN  sphere  drop  calibration  trials.  In 
the  absence  of  ducting,  any  differences  can  generally  be  attributed  to 
plumbing  or  other  isotropic  losses. 


RANRL  has  been  requested  (  ref  3  )  to  produce  programs  suitable 
for  desktop  computers  to  solve  the  inverse  problems  -  (i)  the  calculation 
of  signal-to-noise  required  to  yield  a  given  probability  of  paint  and 
(ii)  the  calculation  and  plotting  of  detection  contours  in  a  multipath 
environment.  The  ensuing  sections  describe  the  algorithms  used  in  the 
programs. 


2.  The  Radar  Equation 


The  power  returned  in  free  space  from  a  target  of  cross-section  a 
is  given  by  the  monostatic  radar  equation 

p  .  PjW*  i. 

r  (4n) »R« 


where  P(  is  the  transmitted  power,  6  is  the  power  gain,  X  is  the  radar 
wavelength  and  R  is  the  target  range.  Multipath,  diffraction  and  other 
environmental  effects  are  accounted  for  by  the  pattern  propagation  factor 
(F)  and  the  atmospheric  loss  factor  (L) 
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The  problem  addressed  in  this  paper  is  the  calculation  of  1.  in  equation 
2,  which  may  be  recast  to  provide  ar,  expression  for  the  mrximur>  sirgle-blrp 
detection  range: 

r  =  Fptc*x’g  10,1>  1-  3. 

,rax  [(4r)JPnD0L 

where  Pt  is  the  peak  transmitted  power  rnd  Pn  if  the  systea.  noise 
power.  Dq  is  the  single-pulse  s igna l-to~noi se  ratio  requiied  to  yield 
the  desired  probability  of  paint  for  a  gi\en  number  of  pulses  integrated, 
false  alarat  rate  and  target  return  statistics.  In  the  absence  of  clutter, 
the  limit  to  signal  detectability  is  governed  by  the  pulse  energy .  so  that 
the  effective  noise  power  P  ,  referred  to  the  artenna,  is  determined  by 
the  transmitted  pulse  width  (x),  the  antenna  noise  temperature  (Ta ) ,  the 
receiving  line  losses  <Lr)  and  the  receiver  noise  figure  (NF) 


pn  “  k^[Ta  +  Tr(Lr-l)  *  LrTc(MF-l)] 


where  Tr  and  Tc  are  the  temperature  of  the  receiving  line  and  290  E 
respectively,  and  k  is  Boltzmann's  constant.  If  the  receiver  noise  band¬ 
width  En  is  used  instead  of  1/x  in  (4)  the  transmitted  power  should  be 
multiplied  by  Bx,  the  time-bandwidth  constant,  to  give  the  effective  S/N 
for  probability  of  detection  calculations.  If  cluttei-to-noise  is  near 
unity,  it  is  convenient  to  assume  that  the  clutter-plus-ncise  variable 
(Pc+Pn)  has  the  same  statistical  distribution  as  receiver  noise,  and  this 
Rayleigh  distributed  total  noise  power  is  used  for  PQ  in  equation  3. 


3.  Required  Signal-to-Ncise 
3,1  Approximate  Formulae 

There  are  numerous  approximate  formulae  in  the  radar  literature 
for  evaluating  paint  probability  from  S/N  (and  vice  versa).  Reasonable 
estimates  of  detection  range  can  be  obtained  using  the  simple  formula 
suggested  by  Neuvyr  (  ref  4  ): 


»o  -  1°  10*  [  fr 


log  PFA 
<log(l/Pd))P 


5. 


where  PFA  is  the  probability  of  false  alarm,  Pd  is  the  probability  of 
paint,  N  is  the  number  of  pulses  incoherently  integrated.  The  detector 
law  is  described  by  the  'constant'  y  which  is  often  given  the  empirical 
value  of  2/3  (  ref  5  )  rather  than  the  asymptotic  limit  of  1/2.  Neuvy 
has  given  heuristic  estimates  of  a  and  f)  for  the  Swerling  and  Marcnm 
(non-fluctuating)  targets  as  shown  in  Table  1. 
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J . 

Case  1 

c  1 

_ 0.. 

1 

1 

Swerling  I  1 

2/3  Cl 

+  2/3  exp(-N/3) } I 

1 

1 

1 

II  1 

1  1 

1/6  + 

exp(-N/3)  1 

1 

III  1 

3/4  [1 

-t  2/3  exp ( -N/ 3  )  ]  1 

2/3 

1 

1 

IV  1 

1  1 

1/6  + 

2/3  exp(- N/3)  1 

Non-fluctuating  1 

1 

-t  2  exp(-N/3)  i 

. 1/  6. 

..  _ i 

Table  1.  Neuvy  parameters  for  Marcum  and  Swerling  targets. 


N 

I 

Swei ling 

III 

Case 

11 

IV 

F.qn  ( 6 ) 

1 

10.5924 

10.3747 

10.5924 

10.3747 

10.5924 

2 

8.1681 

7.9547 

7.9547 

7.9295 

7.9547 

3 

6.8091 

6.5984 

6.5735 

6.5938 

6.7176 

4 

5.8732 

5.6644 

5.6443 

5.6807 

5.6060 

5 

5.1639 

4.9565 

4.9478 

4.9907 

4.8528 

6 

4.5949 

4.3888 

4.3926 

4.4380 

4.2808 

7 

4.1213 

3.9163 

3.9320 

3.9781 

3 .8188 

8 

3.7165 

3.5123 

3.5392 

3.5851 

3.4310 

9 

3.3636 

3.1601 

3 .1972 

3.2425 

3.0968 

10 

3.0511 

2.8483 

2.8947 

2.9391 

2  8030 

20 

1.0674 

0.8689 

0.9763 

1.0105 

0.9560 

30 

-0.03 97 

-0.2360 

-0.0962 

-0.0687 

-0.0798 

40 

-0 . 8042 

-0.99  89 

-0.8384 

-0.8153 

-0.8024 

50 

-1.3861 

-1.5798 

-1.4045 

-1 .3845 

-1  .3574 

60 

-1.8550 

-2.0479 

-1.8613 

-1.8437 

-1.8079 

70 

-2.2471 

-2.4393 

-2.2438 

-2.2280 

-2.1871 

80 

-2.5837 

-2.7754 

-2.5726 

-2.5583 

-2.5143 

90 

-2.8785 

-3.0697 

-2.8608 

-2.8476 

-2.8021 

100 

-3.1405 

-3.3312 

-3.1171 

-3.1049 

-3.0590 

200 

-4.8285 

-5.0169 

-4.7736 

-4.7666 

-4.7370 

300 

-5.7916 

-5.9791 

-5.7225 

-5.7175 

-5.7111 

400 

-6.4664 

-6.6531 

-6.3883 

-6.3844 

-6.3996 

500 

-6.9853 

-7.1715 

-6.9012 

-6.8980 

-6.9324 

600 

-7.4065 

-7.5925 

-7.3180 

-7.3153 

-7.3670 

700 

-7.7611 

-7.9467 

-7.6691 

-7.6667 

-7.7340 

800 

-8.0672 

-8.2526 

-7.9723 

-7.9702 

-8.0516 

Table  2.  Required  signal-to-noise  (dB)  tabulated  for  K-l  to 
800  pulses  integrated  (PFA=0 .000001 ,  P^=0.33). 
Iterative  solutions  (columns  2-5)  used  fitted  data 
in  column  6  as  'first  guess'. 
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These  formulae  are  accurate  to  within  a  few  dE  for  0.1<P,<0.9  and  mod- 

d 

erate  values  of  K.  This  rarge  is  not  adequate  since  (i)  the  95%  contour  is 
often  specified  as  the  required  detectability  contour  and  (ii)  cumulative 
paint  probability  considerations  n.ight  warrant  the  plotting  of  a  Pd  <10*54 
contour.  The  formulae  also  do  net  give  good  agreement  for  1<N<5  which  is 

typical  for  3-D  radar,  nor  are  they  applicable  to  veiy  slowly  fluctuating 
(We instock)  targets.  Expressions  cf  comparable  accuracy  lave  leer  given 
by  AlbersheiO'  (ref  6)  and  Blake  (ref  7)  tor  r.en-f  luctvat  it  g  target... 

3.2  Iterative  Solutions 

The  fcrniulae  described  above  cn  ret  valid  over  a  sufficiently 
large  range  of  P^,  N,  PFA  ar.d  target  scintillation  rate  to  be  used  for 
routine  VCD  calculations,  but  are  sometimes  useful  in  providing  a  starting 
point  for  an  iterative  algorithm.  However,  in  the  unreliable  regions 
(such  as  moderately  large  N  and  F,>90*»  )  numerical  instability  poses  a 
serious  problem.  A  more  robust  starting  point  is  required  which  covers 
the  range  of  radar  fid  target  parameters  likely  to  be  ercour. teied. 

The  method  used  here  is  based  on  the  observation  that  D0,  for  33% 
probability  of  detection  in  gaussian  noise,  is  virtually  independent  of 
the  amplitude  statistics  of  the  target  (  see  figure  1  ).  Regression 
analysis  of  DQ  data  for  Ft,=D.33,  3<N<1000  ,  PFA=10  Swerling  case  II, 
and  non-cohfcrent  integration  yields  the  following  result 

Dq  =  7. 138  +  1.018/ log (N)  -  5.533.1og(N)  6. 

with  D0  in  dB.  Values  for  N=1  and  2  are  evaluated  separately  in  the 
program . 

The  secant  iterative  method  with  equation  6  as  first  guess, 
together  with  the  algorithms  of  reference  2,  were  used  to  produce  the 
data  in  table  2.  Iteration  was  stopped  at  Pd  =  0.33+0.00001.  The 
accvracy  of  equation  6  is  of  the  order  of  the  dependence  on  target 
scintillation  (  10.1  dB  )  at  33%  probability  of  detection.  Results  for 
50%  and  95%  (+  0.001%)  are  shown  in  graphical  form  in  figures  2  and  3. 

4.  The  Radar  VCD 

In  free  space, the  detectability  contour,  or  vertical  coverage 
diagram,  is  determined  by  equation  3  with  F  replaced  by  the  antenna  pattern 
function  f(6) 

where  R0  is  the  detection  rarge  along  boreslght,  and  8  is  the  elevation 
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angle.  This  smoothed  contour  is  also  useful  lor  estimating  mear  cotectitn 
ranges  at  higher  elevation  angles  and  modeiate  sea  states,  in  which  case 
the  multipath  structure  is  washed  out  (see  also  latex  section  on  ship 
motion).  At  lower  elevation  angles,  or  sea  states,  multipath  lobing  must 
be  considered  and  the  detection  center i  becomes 


r-,.„  ■  F  f  o 


t. 


Under  non-ducting  conditions,  the  pattern  propagation  'actor,  F  in  the 
interference  region  is 


F  =  f(ej) ./ 1  +  x2  +  2x  cos  0 


9. 


8nd  may  take  values  between  0  and  2.  The  phase  diffexecce  0  is  the  sum  of 
contributions  fiom  the  geometric  path  difference  between  the  direct  ar.d 
indirect  rays  (  fig  4  ),  and  the  phase  difference  cn  reflect icr  fron  the 
sea  surface  of  the  indixect  ray.  The  reflectivity  paraneter  \  is 


r  p  D  f(02)  10. 

— fWJ 

in  which  D  is  the  divergence  factoi,  i  is  the  roughness  factor,  p  is  the 
dielectric  reflectivity  (the  reflectivity  which  would  apply  if  the  sea 
were  perfectly  smooth )  and  6j  and  ©2  are  *s  shown  in  figure  4. 


For  elevation  angles  near  the  horizon,  and  for  targets  over  the 
horizon,  F  is  calculated  using  diffraction  theory  (or  by  interpolation  as 
described  in  ref  1)  with 

F  =  f(ea) ./U(X) .VTZj) .V(Z2)  11. 

where  X,  Zj,  Z2  are  range,  target  height  and  antenna  height  respectively 
in  natuial  units  and  the  functions  V  and  U  are  gain  functions  described 
in  reference  2. 


4.1  VCD  Envelope 

The  main  features  of  the  VCD  for  Naval  radars  can  be  calcvleted 
using  ray  theory.  The  envelope  of  is  obtained  with  equations  £ 

and  9  with 

0max  =  2nm  (  m=l,2,3 . )  12. 

so  that. 


a 


*Snax  *envel°Pe)  =  f  (0j) .  (l+»)  .Rq 


13. 


-■*-■*-*  v.  '..•via  ^»r  •, t—t 


Within  the  range  of  elevation  ait-  grazing  angles  of  interest, 
the  shape  of  the  VCD  is  thus  dominated  by  f(0)  and  r. 

4.2  The  Roughness  Factor 

The  reflectivity  of  the  indirect  ray  can  be  written  as 

T  =  r.p.e'J®  14. 

where  r  is  the  roughness  factor,  and  p  and  0  are  the  magnitude  and 
phase  respectively  of  the  specular  reflection  coefficient.  Formulae 
for  p  and  0  as  functions  of  grazing  angle,  frequency,  water  temperature 
and  salinity  are  given  in  reference  2  and  are  in  good  agreement  with 
experimental  data.  The  dependence  of  the  roughness  factor  on  grazing 
angle  and  frequency  is  less  straightforward,  and  thera  is  a  paucity  of 
experimental  data. 


Ament  (  ref  8)  has  shown  that  if  the  wave  height  distribution 
is  gaussian,  (variance  a1),  then  the  surface  roughness  will  also  be 


gauss  tan 


where 


,  - 


s  -  2rto .  siny/ X 


This  equation  gives  good  agreement  at  low  grazing  angles  (y), 
but  this  is  to  be  expected  since  r-->l  as  y-- >0.  That  is.  most  models 
will  predict  EjBax~2R0  for  the  lowest  multipath  maximum  over  a  wide  range 
of  frequencies  and  sea  states,  despite  the  lack  of  agreement  for  the  high 
altitude  coverage. 


At  higher  values  of  s  the  reflection  is  net  purely  specular.  The 
additional  diffuse  component  adds  to  the  fluctuation  in  the  pattern 
propagation  factor,  but  net  to  its  average  value.  The  random  component 
of  F  is  not  considered  in  the  program,  but  rather  an  effective  constant 
value  is  assumed.  Reasonable  agreement  for  large  s  is  obtained  using 
the  empirical  expression  given  in  reference  1: 


=  e-l. 2732s 


for  r>0.44 
r<=0.44 


The  program  also  makes  use  of  the  Burling  relationship  between 
significant  wave  height  (F.^^)and  a 

®l/3  =  *  0  * 


'.Vvl  I 


4.3  Antenna  Pattern  Functions 


At  high  elevation  at. ties  and  mctieiatt  sea  states  x--  >0  and  the 
envelope  of  F.^ax  is  dominated  by  the  APF,  as  per  c  qn  7.  Flsewheie,  the 
antenna  pattern  function  of  both  the  direct  and  indirect  rays  a  e 
required  for  the  calculation  of  F.  In  general,  the  APF  needs  to  he 
represented  adequately  out  to  the  first  sidclobe.  The  juc.gran  calculates 
either  (i)  a  cosecant- squared  pattern  or  (  i  i)  a  modified  sin  u/u  far.  bean 
of  the  form 


where  u-d.sin0/A  and  B  is  a  constant  fcr  the  antenna,  which  describes  the 
side-lobe  level  and  aperture  efficiency,  Futfcods  for  calculating  H  from 
the  sidelcbe  level  arc  described  ir  reference  2 


5.  Crapbical  F.epre  sen  t  a  t  ion  of  VCD 

In  the  absence  of  pitch  and  roll  the  VCD  is  independent  of  radar 
azimuth  (neglecting  blind  arcs  and  super st ruct u j e  multipath)  so  that  the 
VCD  can  be  displayed  as  a  2~D  graphical  represser  tat  ion.  Variations  with 
range  and  height  of  paint  probability  or  s  igna 1- to~nc i se  can  be 
described  by  an  arbitrary  numbei  of  grey  scales 

In  figures  5-8,  the  VCD  of  a  F1IF  ail  search  radar  and  a  G~band 
surface  search  radai  are  illustrated  for  two  sea  states.  Gross 
parameters  used  for  the  FHF  radar  are  k  -  0.7  metre,  antenna  height  hj  = 
30  metre,  N-75  pulses  incoherently  integrated,  sea  states  0  and  6,  and  a 
free  space  rarge  of  F.q  =  80  n.  miles  against  a  1  met  re- squa  red  target. 
Parameters  used  for  the  G~band  radar  are  k  0.05  metre,  hj  =  25  metre,  N 
=  5,  sea  states  0  and  3  and  F^  -  17  n. miles.  Standard  atmospheric 
conditions  and  scan-to-scan  (Swelling  case  I  )  target  scintillation  arc 
assumed.  The  grey  scales  correspond  to  paint  probability  regions  P > 9 5% , 
50%<P<95%,  5%<P<50?c  and  P<5%,  and  were  computed  as  described  in  the 
previous  section. 

The  VCD's  were  produced  by  taking  120  cuts  in  height  fcr  200  range 
increments.  Sea  clutter  was  ir eluded  ir  the  noise  calculations,  using 
the  expressions  described  in  reference  2  For  reasons  of  clarity  the 
plots  are  not  truncated  at  minimum  range  and  ran inun  unambiguous  lange  as 
determined  by  the  radiated  pulse  width  and  PRF  respect ir e  ly  .  In  order  to 
resolve  the  structure  in  the  G'band  plots,  calculations  were  carried  out 
only  to  4000  feet  and  limited  to  sea  state  3,  while  the  FHF  calculations 


were  carried  out  to  a  height  sufficient  to  contain  the  5$  probability 
contour . 

The  calculations  for  figures  5-8  took  a  few  minutes  of  processing 
tia>e  on  a  CDC  Cyber  76  mainfrane  computer.  Such  a  program  is  however 
unsuitable  for  small  desktops,  such  as  the  Tektronix  Graphics  System 
specified  in  reference  3,  since  similar  calculations  would  take  2-3  days 
cf  computer  time  for  each  plot.  The  next  section  describes  methods  for 
producing  lire  contour  VCD  s  which  can  be  c;uicklv  computed  on  a  stall 
desktop  machine,  using  a  modification  of  the  progiat  desciibeo  ir  vef  1 

6.  Contour  Plotting 

6.1  Asymptotic  Behaviour  The  Fiat  Faith  Limit 

Fast  algorithms  for  computing  contour  VCD  s  aie  not  easily 
implemented  due  to  the  lack  of  s.ytuictiy  it:  the  multipath  icbiig 
structure.  One  method  often  employed  is  to  pcrfoim  an  approximate 

calculation  of  the  VCD  using  the  flat  ea i th  mul t ipath  lesults  and  to  modify 
the  computed  results  graphically  or  by  scaling  tc  account  foi  cuivatvie 
Lven  this  appioach  howevei  will  not  be  generally  applicable  due  to  the 
additional  geometric  approx  imat ion s  that  rust  be  made 


In  the  flat  earth  limit  D  1  and,  for  horizontally  polarized  UllF 
ndais  at  low  to  modeietc  elevation  angles  r-1  and  the  phase  difference 
on  teflection  is  n  ladians  The  path  difference  between  direct  and 
indirect  rays  for  a  target  at  ground  range  G  is 


(hl‘‘h2)l  +  G* 


(hj-hj ) 2  +  G* 


If  the  free  space  range  of  the  ladar  is  large  compeied  with  the 
sum  of  target  and  antenna  heights  the  path  difference  lor  constructive 
interference  is  given  approximately  as 

6  =  2hlh2  19. 
R  cos8 
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Usually  G>>hj/tar6  fcr  r.aval  search  radars  at  ranges  of  interest . 
so  that  the  VCD’s  produced  from  the  flat  earth  model  are  highly  symmetric 
-  that  is,  only  one  '■'a  leu  lat  ion  of  multipath  geometry  need  be  perftirotd 
for  one  lange  at  each  elevation  anglt  increment,  vith  D"4  scaling  t  <  > 
determine  the  range  for  a  specified  probability  of  detection  and  false 
alarm  rate  (using  the  algorithms  of  section  3).  A  first  order  correct  ion 
for  the  effect  of  the  earth’s  curvature  car  be  included  afttr  the  last 
approximation.  This  is  equivalent  to  assuring  that  the  sea  surface  is 
f’'t  up  to  the  point  of  reflection  so  that  the  final  resilt  requires  c-nlj 
the  transformation  hj— >1^  +  G1/2ae. 

The  algorithms  used  in  the  prograa  described  here  do  net  rely  on  the 
flat  c  r  t  h  model,  however  the  gross  structure  of  the  VCD  can  be  determined 
using  the  procedure  described  above.  Such  a  description  is  therefore  a 
useful  basis  for  a  more  general  algorithm. 

6.2  Spherical  Earth  Model 

-  Dependence  of  F  on  F.ange 

The  flat  earth  model  predicts  that  the  lebe  maxima  for  long  range 
r.aval  radars  occur  (with  F=2  )  at  elevation  angles 

eBax  -  sin*1 (2m-l)X/4h1  m-1,2,3 .  21. 


It  is  clear  from  figures  5~S  that,  in  the  more  general  model,  F-2  at  the  lobe 
maxima  only  at  low  elevation  angles.  This  implies,  not  only  that  twice  the 
fr«e  space  range  will  be  achieved  only  at  lower  altitude,  but  also  that  the 
lobe  spacirg  is  net  uniform.  Equation  21  is,  however,  useful  for  determining 
the  elevation  angle  spacing  required  to  fully  resolve  themultipath  structure. 
In  the  progran  ,  ten  points  are  calculated  per  nominal  lobe  spacirg.  This 
facilitates  both  resolution  of  the  lobes  and  the  ability  to  read  the  radar 
range  from  the  plots  to  a  within  a  few  percent  of  F.  . 

The  behaviour  at  constant  elevation  angle  for  a  spherical  earth 
is  most  easily  demonstrated  graphically.  Figures  9  and  10  are  typical 
plots  of  the  signal  return  for  the  U1IF  radar  deser  ibed  in  the  previous 
section.  The  plots  are  for  a  1  metre-squared  target  at  two  intermediate 
elevation  angles,  separated  by  half  a  nominal  lobe  spacing,  at  sea  states 
0  and  6.  If  the  free  space  range  of  the  radar  is  large  compared  with  the 
clutter  horizon,  algorithms  for  contour  VCD  s  are  likely  to  be  roost 
stable  at  higher  sea  states  due  to  the  reduced  multipath  lobirg.  This  is 
seen  in  the  plot  for  sea  state  6  (figure  10)  where  results  for  both 
elevation  angles  yield  results  which  are  close  to  the  free  space  result. 

Since  the  power  return  decays  as  thi  fourth  power  of  range,  either  a 
scaling  or  iterative  algorithm  should  calculate  the  required 


signal-to-noise  in  one  or  two  iterations  even  il  the  first 
inaccurate . 

The  behaviour  at  low  sea  states  is  less  well  b 
robust  algorithms  are  generally  required.  In  figure  9  t 
at  the  free  space  range  (80  n. miles)  is  qualitatively  di 
two  elevation  angles  selected.  One  curve  corresponds  1 
of  a  lobe  maximum  (at  80  n. miles)  and  R*4  behaviour  is  s 
2.0  times  the  free  space  range.  This  type  cf  behaviour  ;s 
and  implies  that  as  long  as  the  first  guess  for  long 
around  the  free  space  range,  the  detection  range  (S/N=D 
maxima  will  generally  be  easily  computed.  At  very  she 

hopping  occurs  for  constant  elevaticn  angle  calculations  si 
algorithms  may  be  numerically  unstable. 

As  a  corollary,  iterative  algorithms  nay  be  uns 
range  search  radars  such  as  the  G-band  radar  of  the  previc 
the  worst  case  there  may  be  either  no  solution  or  severa 
S/N=T0  at  a  given  elevation  angle.  The  former  case  may  ap 
sea  states  if  the  detection  range  is  of  the  order  of  the  c 
while  the  latter  will  be  worst  at  low  sea  states  when 
lobing  is  most  pronounced. 

6.3  Description  of  Main  Program 

Either  of  two  algorithms  nay  be  selected  in  the  ma 
cases  where  numerical  stability  is  not  a  problem  the  pref 
calculation  is  an  iterative  search  method.  This  algorit 
reliable  results  at  moderate  sea  states  for  naval  radars 
is  large  compared  with  the  clutter  horizon.  This  method 
there  is  a  solution  along  a  selected  elevaticn  angle  and 
is  reasonable  (power  monoton ica 1 ly  decreasing  with  ra 
detection  range  can  be  calculated  with  arbitrary  preci 
sensitivity  cf  C.25  dP  should  be  acceptable  for  most  app 
10  prints  are  plotted  per  nooiinal  lobe  spacing. 

In  order  to  increase  the  probability  of  the  searc 
numerically  stable  region,  the  first  two  corrections  to  tl 
are  scaled  assuming  R“  4  and  P.  ••  decay  law  respectively 
reduce  ringing),  with  subsequent  iterations  independent  of 
I  f  a  solution  has  not  teen  obtained  by  a  specified  number 
the  'solution'  plotted  will  correspond  to  the  minimum  val 
This  method  will  therefore  produce  reliable  ranges  in 
structure  but  not  necessarily  in  the  nulls.  This  loss  of 


of  little  corsequence  (see  figures  5-8)  especially  coo s id* r ing  that  the 
power  in  the  nulls  is  highly  variable  due  to  ship  n.oticn  (discussed  in 
next  section)  ,  wave  height  fluctuations,  atmospheric  inhor ogene i ty ,  and 
other  factors. 

The  second  algorithm  is  a  one-step  scaling  algorithm,  so  that  the 
computed  detectability  range,  at  a  specifitd  elevation  aigle,  will  in 
general  be  reliable  only  if  the  first  guess  is  close  to  the  actual 
detectability  range.  The  scaling  algorithm  in  the  progran  is  optimized 
to  produce  a  reliable  envelope  for  VCD  contour  since  the  first  guess  at 
any  elevation  angle  is  the  range  to  the  lobe  maximum,  given  by  equation  13. 
Since  the  reflectivity  parameter  x  is  a  function  of  R,  it  is  not  known  until 
the  final  solution  is  obtained.  An  estimate  x  is  used,  and  this  provides 
an  estimate 

iax  =  Ro  f<e>.U+*i)  22. 

The  estimate  x^  is  the  value  corresponding  to  the  last  solution 
(x^_j).  With  this  procedure,  I^,ax  will  be  calculable  to  within  a  few 
percent  at  sea  state  0  if  at  least  10  points  are  calculated  per  multipath 
lobe.  At  higher  elevation  angles  and/or  sea  states  the  multipath  lobing 
structure  is  washed  out  (xj — >0)  so  that  this  algorithm  will  be  both 
robust  and  accurate  if  R  is  greater  than  the  clutter  horizon. 

Plots  using  these  algorithms  are  given  in  plots  11-17.  The 
contour  selected  is  C  dB  in  all  cases,  which  corresponds  to  4%  and  52% 
probability  of  paint  for  the  G-band  (N=5)  and  UHF  (N-75)  radars 
respectively  against  a  Swerling  case  I,  1  metre-squared  target 
Figures  11  and  12  used  the  iterative  algorithm  (0.25  dE  sensitivity)  at 
sea  states  0  and  6.  The  contour  VCD's  are  in  good  agreement  with  the 
grey  scale  envelopes  of  figures  5  and  6.  An  additional  plot  for  sea  state 
3  (fig  13)  is  included  to  illustrate  the  gradual,  and  very  significant, 
decrease  in  maximum  height  with  sea  state.  The  sane  parameters  were  used 
to  produce  the  plots  in  figures  14  and  15  with  the  scaling  algorithm. 
Although  this  algorithm  is  optimized  at  the  lobe  maxima,  it  also  gives 
good  agreement  in  the  mid-lobe  region.  The  lower  half  of  each  lobe  is 
generally  well  reproduced  to  much  shorter  ranges  -  this  is  fortuitous 
since  this  is  the  lobe  region  of  interest  for  inbound  air  targets. 

Figure  16  and  17  show  the  results  of  the  scsling  algorithm  for 
the  G-band  contours,  with  the  calculations  stopped  at  an  elevation  angle 
corresponding  to  the  maximum  height  in  the  full  graphical  representations 
(fig  7-8),  Again,  the  0  dB  contour  (4%)  is  consistent  with  the  5%  grey 


scale  envelope  of  figures  7-8.  With  this  radar,  the  pattern  propagation 
factor  has  a  stronger  range-dependi i c r  sc  that  the  iterati\e  algorithm  is 
numerically  less  stable.  The  mar  lesta.'ion  of  this  is  a  VCD  with  a  much 
higher  incidence  of  'bad'  data  prints  (about  5%  of  all  points  calculated). 
At  sea  state  0,  the  lobes  for  this  radar  are  pronounced  but  closely  spaced 
so  that  only  the  envelope  of  the  VCD  contour  is  easily  measured  (siflicient 
reason  for  using  the  quicker  scaling  metluri).  At  hi)  her  sea  stales  the 
effect  of  ship  motion  further  complicates  the  VCD. 

7.  Effect  of  Ship  Motion 

The  motion  of  the  ship  is  most  conveniently  described  in  terms  of 
the  reference  axes  systems  defined  in  figure  18.  A  natural  choice  for 
the  'space-fixed'  axis  system  utilizes  the  mean  sea  surface  as  the  x-y 
plane.  The  ship's  axis  system  is  also  naturally  defined  by  the  effective 
plane  of  symmetry  through  the  keel,  which  is  defined  as  the  x'-z'  plane, 
with  the  y'-axis  passing  through  the  antenna. 

Relative  motion  of  the  two  axes  systems  about  the  vertical  (yaw) 
is  equivalent  to  a  fluctuation  in  the  antenna  rotation  rate  and,  thus, 
the  number  of  pulses  integrated.  Similarly',  the  number  of  hits  per  scan 
is  increased  or  decreased  during  a  rapid  turn.  Although  this  motion  may 
be  of  the  same  order  as  the  normal  antenna  moticn,  the  probability  of 
detection  is  only'  a  weak  function  of  the  number  of  pulses  incoherently 
integrated,  and  can  therefore  be  ignored  in  most  cases.  Motion  of  the 
ship  in  the  x~y  plane  is  also  ignored  since  the  VCD  is  plotted  in  terms 
of  relative  range  which  will  not  change  significantly  during  the  time  on 
target  for  normal  antenna  rotation  rates  and  target  range  rates. 

Pitch  and  roll  can  be  defined  as  the  angles  9  and  0 
respectively  in  figure  18.  The  effects  of  ship  ’rotation*  are  symmetric 
in  pitch  and  roll  unless  a  specific  target  bearing  is  considered.  For  a 
target  on  the  bow,  pitch  has  the  same  effect  as  varying  the  antenna  tilt 
in  the  x,y,z  axis  system  while  preserving  the  polarization  of  the 
radiation.  (The  effect  of  the  vertical  motion  of  the  antenna  during 
pitch/roll  is  discussed  later  in  the  section  on  heave).  Pitch  or  roll 
can  be  significant  compared  with  vertical  beamwidth,  even  for  wide 
beamwidtb  search  radars.  This  not  only  has  the  effect  of  increasing  the 
maximum  angle  of  the  main  antenna  lobe  VCD,  but  also  gives  rise  to 
calculable  fluctuations  in  the  pattern  propagation  factor  at  moderate 
elevation  angles  due  to  variations  in  the  APF  of  the  direct  and  indirect 
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rays.  This  is  shown  in  figure  19,  where  the  antenna  tilt  was  allowed  to 
vary  randomly  between  zero  and  10%  of  the  vertical  beamwidth. 

The  effect  of  roll  tor  a  target  along  zero  relative  bearing  also 
has  a  small  a  ■/.  imu  lb  f  luct  ua  t  i  1 1.  effect  for  targets  at  moderate  elevation 
angles.  (This  effect  is  not  normally  significant,  or  else  it  would 
provide  an  elegant  method  of  de  t  erm  ir.  ing  target  elevation  using  a  2~D 
radar.)  A  second  effect  for  this  relative  geometry  is  that  polarization 
of  the  radiation  in  the  space-fixed  system  is  partially  converted  to  the 
opposite  sense.  The  detecting  antenna  is  only  concerned  with  the  polariz¬ 
ation  in  the  ship's  axis  system,  so  that  this  is  only  manifested  through 
multipath  effects.  Six  degrees  of  roll  converts  only  1%  of  the  radiation 
to  the  opposite  polarization.  At  grazing  incidence  for  the  indirect  ray, 
the  phase  difference  anc  reflectivity  arc  near  rt  and  1  respectively  for 
both  polarizations,  so  that  the  lowest  lobe  is  virtually  independent  of 
polarization.  At  moderate  elevation,  the  magnitude  and  phase  for  vertical 
polarization  may  be  sufficiently  different  to  put  maxima  at  the  elevation 
angle  of  a  minimum  for  the  opposite  polarization.  Since  the  detection 
range  is  of  the  order  of  (S/K'°'is,  the  nulls  for  a  1%  polarization 
change  may  be  filled  to  about  a  third  of  the  range  for  the  adjacent  lobe 
maxima . 

In  the  case  of  heave,  the  effect  is  simply  related  to  the  ratio 
of  heave  to  mean  antenna  height  above  sea  level.  A  simulation  of  the 
effect  of  heave  is  shown  in  figure  20  (one  run  only), where  the  antenna 
height  was  uniformly  distributee1,  over  the  nominal  heave  dimension.  The 
effect  on  the  lower  lobes  is  only  slight,  but  heave  has  the  effect  of 
filling  in  the  upper  lobes  and  thus  reducing  the  mean  detection  range 
along  a  lobe  maximum.  If  the  mean  of  many  simulation  runs  is  calctlated 
for  each  elevation  angle,  the  nulls  at  the  n-th  lobe  will  be  completely 
filled  if  n  is  of  order  (antenna  height/2. heave)  or  greater,  decreasing 
in  effect  with  decreasing  elevation  angle. 

Ship  heave  will  also  affect  the  probability  of  detection  by  its 
effect  on  the  target  fluctuation  statistics.  In  the  case  of  slowly 
fluctuating  (Weinstock)  targets,  hc;>ve  will  modify  the  scintillation  to 
scan-to-scan  (Swerl ing  case  I)  at  higher  elevation  angles  but  will  have  a 
reduced  effect  on  the  amplitude  statistics  at  lower  atgles.  Targets  with 
Swerling  case  I-IV  statistics  will  not  be  affected  to  the  same  extent 
since  ship  motion  is  negligible  on  a  pul se-to-pulse  timescale  for  typical 
PRF’s. 


8.  Frequency  Ag  i  I  i  t  \  and  Diversity 

The  radars  discussed  ai"  \  e  are  assumed  to  have  a  t  ran  sir.  it  t  id 
frequency  bandwidth  of  crdei  il  the  hi'ii  st  cf  the  pulsewidth  -  typicaliy 
1  MHz.  With  pulse  compression  :u  -la  .<  it  is  the  inverse  of  the  compressed 
pulsewidtb.  In  addition,  the  conic  *  j  equercy  may  be  tuned  over  a  range 
of  several  percent  -  typically  1 1  r,  s  o:  h‘Jl  z .  The  single-pulse  bandwidth 
is  small  compared  with  the  centre  frequency,  and.  so  tas  nc  significant 
affect  on  the  VCD,  while  the  tunability  of  the  radar  set  simplj  changes 
the  number  of  lobes  that  fit  into  the  ancenna  main  lobe  at  fairly  long 
time  intervals.  (The  frequency  term  tu  the  radar  equation  can  be  assumed 
to  be  a  constant,  sir.ct  the  detection  range  varies  as  the  square  root  of 
the  wavelength  . ) 

Frequency  agility  has  an  analogous  effect  on  the  radai  VCD.  The 
detectability,  averaged  over  all  transmitted  frequencies,  may  be  the  same 
as  a  simple  tunable  radar,  however  on  a  si  an- to-scan  timescale  the  radar 
VCD's  are  quite  different.  The  elevation  angle  to  the  n-th  lobe  is 
approximately  proportional  to  the  ratio  of  the  transmitted  wavelength  to 
the  antenna  height.  If  the  agility  was  random,  and  on  a  scan-to-scan 
basis,  the  effect  on  the  VCD  would  be  similar  to  the  heave  simulation  in 
figure  20.  That  is,  the  lower  lobes  would  be  unaffected  but  the  power  in 
the  upper  lobes  would  be  fluctuating  about  the  free  space  level,  Pulse- 
to-pulse  random  agility  yields  the  same  mean  detection  range  but  the 
fluctuations  are  averaged  out  in  the  integration  process.  A  second 
effect  of  random  pul se-to-pul se  frequency  agility  is  on  the  number  of 
independently  fading  signal  groups  per  scan,  or  alternately  the  number  cf 
degrees  of  freedom  of  the  equivalent  Chi  square  target.  If  F  frequencies 
are  transmitted  per  scan  with  sufficient  separation  to  have  independent 
echoes,  then  a  target  which  is  represented  as  having  2K  degrees  of 
freedom  for  the  fixed  frequency  radai  has  up  to  2Kh  degrees  of  freedom 
for  the  pulse-to-pul se  frequency  agile  radar  (  K=F,N,2F  and  2N  for 
Swerling  cases  I, II  III  and  IV  respectively  ). 
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ANNEX  A 

Calculation  of  Detection  Threshold 
For  Fixed  Threshold  Detectors 


ft 

ft 

ft 
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1201  REM  imHiHHiimHHmimmHommoHimHmmHHiHH 

1205  REM  UTILITY  ROUTINE  TO  CALC  SUM  N0=0  TO  N-1  OF  Y0*M.e“-Y0/K ! 

1206  DEF  FN  LGT ( X ) =L0G ( X ) /LOG (10) 

1210  IF  N>1  GOTO  1215 

1211  Y2=EXP(-Y0) 

1212  RATIO=1.0 

1213  GOTO  1280 

1215  NO =0: ANSWER =0 

1216  LIMIT=1.0E35 

1217  IF  Y0>1  THEN  Y1 =Y0»1C“(37-  FN  LGT ( YO )): ELSE  Y 1  - 1/Y0*10“ (37+  FN  LGT(YO) ) 
1220  FACT0R=-L0G(Y1) 

1225  REM  START  OF  MAIN  LOOP 
1230  FACT0R=FACT0R+L0G(Y1 ) *2 
1235  Y 1 = 1 /Y 1 

1240  IF  N0=0  THEN  Y2=Y1:ELSE  Y2=0 
1245  REM  START  OF  INNER  LOOP 
1250  N0=N0+1 

1255  Y1=Y0/N0*Y1 

1260  Y2=Y2+Y1 

1265  IF(N0<(N-1))  AND(Y2<LIMIT)  GOTO  1245 

1 270  ANSWER= ANSWER+EXP ( LOG ( Y2 ) +FACT0R- YO ) 

1275  IF  N0<(N- 1 )  GOTO  1225 

1276  IF  Y1>0  THEN  RATIO=EXP(LOG(ANSWER)+YO-FACTOR-LOG(Y1 ) ) 

1278  Y2=ANSWER 

1280  RETURN 
1290  REM 

1320  REM  START  OF  MAIN  ROUTINE  TO  CALCULATE  THRESHOLD  (YO)  FROM  N  AND  PFA 
1330  Y0=-1«L0G(PFA) 

1335  RATI0=1.0 

1340  IF  N<=1  THEN  1410 

1350  REM  First  estimate  for  YO 

1360  YO  =  (SQR ( -LOG ( PF A ) ) +SQR ( N ) - 1 ) »SQR ( -LOG ( PFA ) ) +N-SQR ( N ) 

1370  GOSUB  1201 

1380  D0=L0G(Y2/PFA)*RATI0 

1390  Y0=Y0+D0 

1400  IF  ABS(DO/YO)=>3.0E-7  THEN  1370 

1410  REM  YO  IS  THRESHOLD  FOR  FIXED  THRESHOLD  DETECTOR 

1420  RETURN 

1425  REM 


Calculation  of  Pd  for  Marcum,  Swerling  and  Weinstock  Targets 


1425  REM 

1430  REM  CALCULATE  Pd  FROM  S/M  (X  dB) , CHI-SQUARE  PARAMETER  (K),N  AMD  YO 
1440  ELIMIT=1/(10*N):Y2=PFA 

1443  Y1=PFA/RATI0 

1444  N0=N 

1450  X=10*(X/10) 

1460  S1=Y2 

1470  X1=(K/(K+N«X))'K 

1480  X2=X 1 

1490  X3=X*N/(K+N*X) 

1500  D1=K-1 

1510  M=1 

1520  R1 =0 

1525  IF  X2=0  THEN  GOSUB  1661 

1526  REM  RESUME  AFTER  UNDERFLOW  LIMIT  REACHED 
1530  L 1 =R 1 

1540  Y1=Y1/N0*Y0 

1550  Y2=Y2+Y1 

1560  S 1=S1+Y 1*( 1-X2) 

1570  El =( 1-X2)*( 1-Y2) 

1580  R1=S1+E1 

1590  X1=(D1+M)/M*X1*X3 

1600  X2=X2+X 1 

1605  IF(X2>1 )  THEN  X2=1 

1610  M=M+1 

1620  N0=N0+1 

1630  IF  ABS (1-L1/R1)=>3. 0E-7  THEN  1530 

1640  IF  E1=>ELIMIT  THEN  1530 

1650  REM  R1  IS  THE  PROBABILITY  OF  DETECTION 

1655  RETURN 

1656  REM 

1660  REM  SUBROUTINE  TO  SCALE  FOP  UNDERFLOW 

1661  X4  =0 : X5  =K« ( LOG ( K ) -LOG (K+N*X) ) 

1662  X6=LOG(X3) 

1663  REM  JUMP  HERE  TILL  LIMIT 

1664  Y1=Y1/N0*Y0 

1665  Y2=Y2+Y1 

1666  S1=S1+Y 1 

1667  X4=X4+X6+LOG(D1+M) -LOG(M) 

1668  X1=EXP(X5+X4) 

1669  X2=X2+X1 

1670  M=M+1 

1671  N0=N0+1 

1672  IF  X2=0  THEN  GOTO  1663 

1673  RETURN 
1675  REH 


ANNEX  C 


Calculation  of  Required  Signal-to-Noise  (Do) 


8000  REM  mmHmiHHmmmHiimimmmmHHHmgimiHt 

8001  REM  routine  to  solve  roots  of  f(x)=R1  by  secant,  method 

8002  REM  SNdBn  IS  CURRENT  ESTIMATE  IN  dB  OF  S/N  REQUIRED  FOR  Pd  =  PROB 

8003  DEF  FK  LGT(X) =L0G(X)/L0G( 10) 

800H  LN=  FN  LGT(N) 

8005  REM 

8006  PROB =0.95 

8007  PROBLIMIT=0. 00001 

8008  REM  ITERATION  WILL  BE  STOPPED  WHEN  SNaB3  =  PROB  +/-  PR0BL1MIT 

8010  REM  FIRST  ESTIMATE  OF  REQUIRED  S/N  IS  FIT  OF  SOLUTIONS  FOR  PD=0.33 

8011  IF  N<3  THEN  SNdB3=10. 5924-8. 74904*LN 

8012  IF  N>=3  THEN  SNdB3=7. 138+1 .018/LN-5.353*LN 

8013  SNdB3=SNdB3+4.343*  FN  LGT(PFA/1 .0E-6)/  FN  LGT(PFA) 

80U  REM  SECANT  METHOD  SEEDED  WITH  2  POINTS  STRADDLING  Pd=0.33 
8015  SN'dBI  =SNdB3-1 .0  :X=SNdB1  :G0SUB  1430 :PROB1 =R1 

8020  SNdB2=SNdB3+1.0:X=SNdB2:G0SUB  1430:PR0B2=R1 

8025  REM 

8028  IF(PROF2=0)  AND(PR0B1=0)  THEN  SNdB3=SNdB3+0 . 25 :GOTO  8015 

8029  IF(PR0B2=1)  AND(PR0B1=1)  THEN  SNdB3=SNdB3 -0.25: GOTO  8015 

8030  SL0PE= ( SNdB2-SNdB 1 ) / ( PROB2-PROB 1 ) 

8035  TESTSL0PE=(PR0B-PR0B2)*SL0PE 

8036  IF  TESTSLOPE>3  THEN  SNdB3=SNdB2+  FN  LGT(TESTSLOPE) rGOTO  8050 

8037  IF  TESTSLOPE<~3  THEN  SNdB3=SNdB2-  FN  LGT(-TESTSLOPE) :GOTO  8050 

8040  SNdB3=SNdB2+TESTSL0PE 

8050  X=SNdB3 :GOSUB  1430:PROB3=R1 

8055  REM  SNdB3  IS  CURRENT  ESTIMATE  OF  Do  FOR  Pc=100«PROB3  % 

8060  IF  ABS ( PR0B3-PR0B ) <PROBLIMIT  THEN  GOTO  8100 

8070  SNdBI =SNdB2 :PR0B1 =PR0B2 

8080  SNdB2=SNdB3:PR0B2=PR0B3 

8090  GOTO  8030 

8100  PRINT  "ITERATION  STOPPED  AT  ";SNdB3;"  dB  =  " ; 100*PROB3 ; " 

8110  RETURN 

9000  REM 
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ANNEX  D 


Calculation  of  Paint  Probability  for  Non-fluctuating  Targets 


REM  «*«*  PROBABLILTY  OF  DETECTION  FOR  KON-FLUCTUATING  TARGETS*** 
REM  CALCUATION  OF  Pd  AT  :  Sigr.al-to-Nci se  =  X  OB 


REM  Probability 

REM  N  Pulses  Non 

REM  Fixed  Thresh 

REM 

Y2=PFA 

Y1=PFA/RATI0 

N0=N 

X=10"(X/10) 

S 1  =Y2 

X3=N*X 

X1=EXP(-X3) 

X2=X1 

X6=LOG(X3) 

M=1 
R1  =0 

REM  TEST  FOR  UNDERFLOW  CONDITION 
IF  X2=0  THEN  GOSUB  1661 
L 1  =R1 

Y1=Y1/N0*Y0 
Y2=Y2+Y 1 
S1=S1+Y1*( 1-X2) 

El =( 1-X2) *( 1-Y2) 

R1=S1+E1 
X1=X3/M*X1 
X2=X2+X 1 
M=M+ 1 
N0=N0+1 

IF  ABS( 1 -L1/R1 ) =>3 .OE-7  THEN  1530 

IF  E1=>0.001  THEN  1530 

REM  R1  IS  THE  PROBABILITY  OF  DETECTION 

RETURN 


Probability  of  False  Alarm  =  PFA 
N  Pulses  Non-coherently  Integrated 
Fixed  Threshold  =  Y0 
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A  N'fvTX  E 

Main  Program  for  Contour 
VCD  calculations. 

REM  DEFINE  FUNCTIONS  REQUIRED  FOR  GEOMETRY  ROUTINES 
DEF  FN  LGT(X) =0.43429448*LCG(X) 

DEF  FN  ASN(X) =ATN(X/SQR( 1-X*2) ) 

DEF  FN  SL (X ) =SQR( (HH2-HK1 )'2+4* ( A1+HK1 ) *( A1+HIP) *(SIN(X/(2*A1 ) ) ) *2) 

DEF  FN  GR(X) =2*A1*  FN  ASN(SQR( ( (X*2- (BH2-HH1 ) “ ?)/ ( 4*( A1+HH1 ) * ( A 1 +HH2) ) ) 
DEF  FN  EL (X ) =  FN  ASN( (2* A 1 « (HH2-HH1 ) +HH2* 2-HH T 2-X*2) / ( 2* ( A 1+HH 1 )«X) ) 
DEF  FN  IND(X) =  FN  ASN( (2*A 1 *HH 1 +HH l“2+X~2) / (2* ( A1 +HH1 ) *X) ) 

REM  VCL  ALGORITHM  IS  ITERATIVE  OR  R*4  SCALING  (  I  OR  S  ) 

REM  30  MULTIPATH  LOBES  CALCULATED  USING  DEFAULT  PARAMETERS  (POINTS* =300 


REM  START  OF  MAIN  PROGRAM  j 

R7=SQR(2*A1/M1 )*(SQR(H1 ) ) 

MiINELEV  =- 1 . 0*  FN  ASN(H1M/R7+R7/2/A1 ) : MAXELEV=TILT+ 1 . 1  *B  1 
PRINT"MIN  ELEVATION  =  " ;57 . 3*MINELEV DEGREES" 

DBL I MI T  =0 . 25 : REM  0.25  dB  RESOLUTION  FOR  ITERATIVE  ALGORITHM 
PRINT"CALCS  CARRIED  OUT  TO  " ;57 . 3*MAXELEV DEGREES" 

INCELEV=0 . 1*  FN  ASN(W/(2*K1 ) ) 

ELTHETA( 1 ) =0 . 70*KINELEV 

PRINT"USING  LOOKUP  TABLE  PROVIDED, WHAT  IS  THE  SIGNAL-TO -NOISE  RATIO" 
PRINT "REQUIRED  FOR  THIS  RADAR'S  VCD  CONTOUR" : INPUT  THRESH 
LINP1 =10* (PI/ 10) 

QTHRESHslO*  (-THRESH/40)  :R0=R0*QTH R’ESH 
FOR  1=1  TO  POINTS* 

FIRSTSLANT=ABS(R0*F1*( 1+R1*F2/F1) ) 

IF  1=1  THEN  GOTO  2220 

ELTHETA ( I ) =ELTHETA ( I- 1 ) +INCELEV 

IF  ELTHETA ( I ) >MAXELEV  THEN  I=POINTS*  :GOTO  2400 

THETA=ELTHETA(I) :GOSUB  920 

REM  First  estimate  of  range 

SLANT=FIRSTSLANT 

PRINT"  FIRST  SLANT  =  " ; SLANT 

LOWLIMIT= 10 

ITERATIONS 

REM  ark**  Compute  first  estimate  of  altitude 
IF  ALG$="I"  GOTO  2255 
IF  ITERATI0N=2  GOTO  2375 
HH1=H1M 

H9(I) =SLANT~2+2*SLANT*( A1+HH1  )*SIN(  ELTHETA  (I) )  +  (A1  +HH1 )  “2 
H9(I)=M1*(SQR(H9(I))-A1) 

PRINT"HEIGHT  =  ";H9(I);"  FEET" 

HH 1=H 1M:HH2=H9 ( I ) /Ml :G(I)=  FN  GR(SLANT) 

H2=H9(I) 

G8=G(I) 

ITERATION: ITERATION+1 
R7=SQR(2«A1/M1)*(SQP(H1)+SQR(H9(I))) 

REM  Compute  target  multipath  geometry 
GOSUB  2680 

REM  Compute  clutter  return  for  itb  point  -  C7(I) 

GOSUB  3400 

REM  Calculatiori/interpolation  of  pattern  propagation  factor  -  F(I) 
GOSUB  2470 

REM  radar  equation 

F(I) =K6+F(I) +RCS-40*  RN  LGT (SLANT) -2*L3*SLANT 


J .  »'!  •'  J.  V' 


•  1 


21 


Fain  program  (cor.t'd) 


2352  NOISE=10*  FN  LGT (LTNP 1+( 10* (C? ( I ) / 10) ) ) 

2353  INCSL  AKT=SL  AM-OLDSL  AM 

235^4  OLDEXCESS=EXCESS:OLDSLANT=SLANT 

2355  EXC ESS  =  F ( 1 ) -HOI £E : I F  ABS(EXCESS-THRESH) CDBLI FIT  THEN  GOTO  2375 

2356  I NC EXCESS =EXCESS-OLDEXCESS 

2357  IF ( ABS ( EXC ESS-THRESB ) > ABS ( LOlrfLlKIT ) )  THEN  GOTO  2359 

2358  LOWLIMI T  =  EXCESS-TH  RESH : LOWG  =G ( I ) : LOWR=OLDSLANT : L0WH=H9 ( I ) :LOWP=F(I) 

2359  IF  ITERATIONC 1G  THEN  GOTO  2365 

2360  EXCESS =LOWLIMIT+TRRESH : G ( I ) =L0WG : OLDSLANT=LOWP : H9 ( I ) =LOWH : F ( I ) =LCWP 

2361  GOTO  2375: REM  END  ITERATION  AND  USE  SMALLEST  EXCESS  IN  VCD 

2362  IF  ITERATIONS  THEN  SLANT=SLANT*QTHRESH* ( 10*  (EXCESS MO) ) 

2363  IF  ITERATIONS  THEN  GOTO  2369 

2365  IF  ITERATIONS  THEN  SLANT=SLANT* (QTRPESH* ( 10' (EXCESS/40) ) )*0.5 

2366  IF  ITERATIONS  THEN  GOTO  2369 

2367  SLANT=SLANT+( THRESH-EXCESS ) *INCSLANT / I NCEXCESS 

2368  IF  SLANT<0  THEN  SLANT=2*R0*RND( 1 )+0 .081 1*W5 

2369  PRINT"SLANTs  " jOLDSLANT ELEV=  " ;57.3*ELTHETA(I) FXCESS=  EXCESS; 

2370  PRINT: GOTO  2250 

2375  RSLANT(I) =OLDSLANT 

2380  REM  Printout 

2390  GOSUB  2840 

2395  IMAXrl 

2400  NEXT  I 

2450  REM 
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Figure  2.  Slgnal-to-noise  required  for  50X  paint  probability 
PFA  -  0. 000001,  Sterling  case  E— IV  targets. 
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Figure  3.  Signal-to-noise  required  for  95X  paint  probe 
PFA  c  0.000001,  Sweriing  case  I-IV  targets. 


Figure  6.  Vortical  coverage  diagram  for  UHF  radar. 
Sea  state  «  6 

Number  of  pulses  integrated  :  75 
Target  j  1  m‘2.  Sterling  case  I 
Grey  scales  i  4 
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Figura  10.  Signal  returns  for  UHF  radar  at  constant  elevation  angle 
of  7.0  end  7.35  degrees.  Sea  state  6. 
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Figure  11.  Contour  VCD  for  UHF  rodar. 

Sea  state  :  0 

Signal-to-noise  *  0  dB  (+  0. 25) 

Algorithm  i  Iterative 
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Figure  14.  Contour  VCD  for  UHF  radar. 
Soa  state  :  0 
Signal-to-noise  i  0  dB 
Algorithm  i  Scaling 
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Figure  16.  Contour  VCD  for  G-band  radar  (  low  angle  coverage  only  ). 
Sea  stGte  :  0 
Signal-to-noise  t  C  dB 
Algorithm  :  Scaling 
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Contour  VCD  for  G-band  radar  Clour  angle 
Sea  state  :  3 
Signal-to-noise  t  0  dB 
Algorithm  :  Scaling 
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Figure  19.  Simulation  of  the  effect  of  varying  antenna  tilt  for  UHF  radar 
Sea  state  t  3 
Signal-to-noise  *  0  cffl 

Antenna  tilt  :  Random  between  -10X  and  +  10X  of  vert  beamwidth 
Algorithm  :  Scaling 
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Figure  20.  Simulation  of  effect  of  varying  antenna  height  for  UHF  radar 
Sea  state  >  3 
Signal-to-noise  i  0  dB 

Antenna  height  fluctuation  »  ±  2.52  of  nominal  height 
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Figure  9.  Signal  returns  far  UHF  radar  at  constant  elevation  angle 
of  7.0  and  7.35  degrees.  Sea  state  0. 
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